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' Abstract 

o ■ 

' Small scale turbulent dynamo action in compressible transonic turbulence is discussed. It 

is shown that the critical value of the magnetic Reynolds number displays a bimodal behavior 
and changes from a typical value of S5 for small Mach numbers to about 80 for larger Mach 
numbers. The transition between the two regimes is relatively sharp. The direct simulations 
are then compared with simulations where shocks are captured using a shock viscosity that 
becomes large at locations where there are shocks. In the presence of shear it is shown that 
large scale dynamo action is possible. 
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1 Introduction 
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There are two main classes of dynamos, small scale and large scale dynamos. This distinction 
^ is not completely unambiguous because in the presence of shear any small scale field will 
*0 ■ attain a large scale component in the direction of the shear. Therefore, the class of small 
scale dynamos is here confined to the case of nonhelical, isotropic and homogeneous turbulent 
^ . flows (32). With shear, however, even nonhelical dynamos can produce large scale fields O. 
Oh| Members of the large scale dynamo class include turbulent flows with sufficient amounts 

Q . of helicity and/or shear such that a large scale held is generated. Here, a large scale field 
i—{ is the component that survives after averaging. The averaging procedure has to be defined 
c/2 1 appropriately and depends on circumstances and also on what kind of field is generated. In 
the presence of a shear fiow, Uq{x, z) = (0, [/q, 0), for example, a useful average (denoted by 
an overbar) will be 



B^^l Bdy. (1) 



H ' Jo 
CZ . . . . . 

■ " This definition will obviously preclude the study of nonaxisymmetric fields. In the case of 

helical turbulence, we have essentially an a effect, i.e. the turbulent electromotive force has a 
field-aligned component and u x b ■ B ^ Q. Here, b = B B is the fluctuating magnetic field 
and u = U — U is the fiuctuating velocity. In unbounded space, for example, the mean field 
can have any orientation. Even in a triply-periodic domain there are still three ultimate field 
configurations, corresponding to Beltrami waves with variation in any of the three coordinate 
directions. The appropriate mean field is then best defined as a two-dimensional average over 
the other two coordinate directions. 

Ideally, when defining averages we want to conform with the Reynolds rules (e.g. Krause 
& Radler 1980). In particular, we want to make sure that the average of an average gives 
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the same average (which is not the case for a running mean), and that average and derivative 
operators commute (not the case for averages over non-periodic directions, including time 
averages). We also want the average of a product of fluctuating and mean quantities to 
vanish (which is not the case for spectral filtering). Therefore, for many practical purposes, 
equation Q is the preferred choice, avoiding any of the aforementioned problems. 

A more formal characteristics of a successful large scale dynamo is therefore one where 
the ratio 

q^{B')/{B')»R-,\ (2) 

Here, angular brackets denote volume averages. Values of q around 0.7 are typical under ideal 
conditions (see Spiral galaxies tend to have q > 0.2 (e.g. Beck et al. 1996). For the sun 
the value of q is unclear, but we would still classify it as a large scale dynamo even if q was 
as small as 0.01, say. 

In the following we give a brief overview of recent progress in the fields of small scale and 
large scale dynamos. In this review we focus on numerical results. 

2 Small scale dynamos 

Small scale dynamos are generally much harder to excite than large scale dynamos. In fact, for 
unit magnetic Prandtl numbers (PrM =v/rj— 1) the critical value of the magnetic Reynolds 
number, i?ni = u^ms/iilkf), is 35 for the small dynamo (e.g., Haugen et al. 2004a) and only 
1.2 for fully helical dynamos (Brandenburg 2001). Here we should emphasize that there is an 
advantage in defining i?ni with respect to the wavenumber instead of the forcing scale (which 
would make R^a larger by 27r) or even the scale of the box (which would make Rm larger 
by the scale separation ratio). The advantage is that with our definition (which is actually 
quite common in the forced turbulence community) the value of i?,n can be regarded as a 
reasonable approximation to rjt/r], where rjt is the turbulent (effective) magnetic diffusivity 
and 77 is the microscopic value. 

In the following we summarize what is now known about the energy spectra in the linear 
and nonlinear regimes and what happens in the presence of shocks. 

2.1 Kazantsev and Kolmogorov spectra 

In the kinematic regime, the magnetic field is still weak and so the velocity field is like in the 
nonmagnetic case, with the usual Kolmogorov spectrum followed by a dissipative subrange. 
Due to an extended bottleneck effect^, only marginally indicated; see Fig. ^ for a simulation 
by Haugen et al. (2004a) where PrM = 1 and i?m — 600. During the kinematic stage, the 
magnetic field shows a clear fc"^/^ spectrum that is characteristic of the Kazantsev (1968) 
spectrum that was originally only expected in the large magnetic Prandtl number limit, i.e. 
for Privi = v/r] ^ 1. During the kinematic phase the spectral magnetic energy grows at all 
wavenumbers exponentially in time and the spectrum remains shape-invariant. 

As the magnetic energy increases, the spectrum arranges itself underneath an envelope 
given by the original kinetic energy spectrum. During this process, the kinetic energy de- 
creases by a certain amount and, above a certain wavenumber, the field can be in super- 
equipartition with the velocity; see the right hand panel of Fig. ^ for ^ high resolution run 

'^For details regarding the bottleneck effect see the recent paper by Dobler et al. (2003), where the difference 
between the fully three-dimensional spectra and the one-dimensional spectra available from wind tunnel 
experiments is explained. 
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Figure 1: Left: early spectra of kinetic and magnetic energy, normalized by ^u^^^/ki, during 
the kinematic stage of run D2. [Adapted from Haugen et al. (2004a).] Right: magnetic, kinetic 
and total energy spectra. 1024'^ meshpoints. The Reynolds number is Urms/(j^fcf) ~ 960. 
[Adapted from Haugen et al. (2003).] 
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Figure 2: Magnetic energy spectra for runs with magnetic Prandtl numbers ranging from 
0.3 to 30. [Adapted from Haugen et al. (2004a).] 



with 1024"^ meshpoints (Haugen et al. 2003). These spectra are, as usual, integrated over 
shells in k space and normalized such that / EKdk — and / Eudk = ^{B^)/^o- The 

magnetic energy displays a nearly flat spectrum in the range 1 < /c < 5, peaks at A: « 5, and 
begins to exhibit an inertial range in 8 < A; < 25, followed by a dissipative subrange over one 
decade. In the inertial range Euik) / EK{k) is about 2.5. 

At larger magnetic Prandtl numbers one begins to see the possible emergence of a k~^ 
tail in the magnetic energy spectrum; see Fig. |21 The k~^ tail has recently been found 
in large magnetic Prandtl number simulations with an imposed magnetic field (Cho et al. 
2002). The k~^ spectrum has its roots in early work by Batchelor (1959) for a passive 
scalar and Moffatt (1963) for the magnetic case. In the low magnetic Prandtl number case, 
evidence is accumulating that the critical magnetic Reynolds number continues to become 
larger (Schekochihin et al. 2004, Boldyrev & Cattaneo 2004, Haugen et al. 2004). 
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Figure 3: Critical magnetic Reynolds number ReM.crit as a function of Ma for simulations 
with PrM = 1 (left) and PrM = 5 (right). Note that RcM.crit depends strongly on Mach 
number for Ma « 1. The simulations with shock-capturing viscosity give approximately the 
correct growth rates. The simulations that provide these data points have resolutions ranging 
from 64'^ to 512'^ mesh points. [Adapted from Haugen et al. (2004b).] 
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Figure 4: Grey (or color) scale representation of V • ix (left) and (right) in an xy cross- 
section through z = for Ma =1.1 using constant viscosity (Run 3a of Haugen et al. 2004b 
with 512^ meshpoints). [Adapted from Haugen et al. (2004b).] 



2.2 Do shocks kill dynamos? 

In the interstellar medium the flows are generally supersonic and in places even highly super- 
sonic. Naively, one could think of highly supersonic turbulence as being nearly irrotational, 
in which case the dynamo should actually become much more efficient and the growth rate 
should increase with Mach number (Kazantsev et al. 1985). This is now known not to be 
the case: in a mixture of irrotational (a; = V x m = 0) and solenoidal (V ■ u = Q) flows the 
critical values of the magnetic Reynolds number for the onset of dynamo action does actu- 
ally increase as a function of the ratio <t = ((V • uY) / (a;^) (Rogachevskii & Kleeorin 1997). 
However, the ratio a is found to stay finite in the limit of large Mach numbers (Porter et al. 
1998, Padoan & Nordlund 1999). This explains the recent finding that minimum magnetic 
Reynolds number for dynamo action displays a bimodal behavior as a function of the Mach 
number Ma = Urms/cs (Haugen et al. 2004b). In fact, for PrM = 1, they find i?m.ciit ~ 35 for 
Ma ^ 1 and i?m,crit ~ 80 for Ma ^ 1. For PrM — 5, the critical values are a bit lower (25 
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Figure 5: Energy spectra for Runs la and Ic of Haugen et al. (2004b). The dotted lines give 
the result using shock-capturing viscosity. [Adapted from Haugen et al. (2004b).] 

and 50 respectively); see Fig. 13 Indeed, in these simulations the ratio a is always about 1/4. 
Visualizations of the magnetic field show that the effects of shocks is rather weak; see Fig. 0] 
where we compare both cases. 

It is worth noting that the use of a shock capturing viscosity {not used in the calculations 
presented in Fig. ^ seems to reproduce the critical values of the magnetic Reynolds number 
rather well; see the left hand panel of Fig.O In Fig. Owe compare kinetic and magnetic energy 
spectra for the constant and shock-capturing viscosity solutions. The two show excellent 
agreement at low wavenumbers. At high wavenumbers there are major differences. The direct 
simulations show an extended diffusion range which is rather from the diffusive subrange 
in incompressible simulations (e.g. Kaneda et al. 2003). The large extent of this diffusive 
subrange in supersonic turbulence is obviously the reason for the tremendously high resolution 
needed in the direct simulations. Fortunately, as far as dynamo action is concerned, this 
extended subrange does not need to be fully resolved and it suffices to just cut it off with a 
shock-capturing viscosity. 

3 Large scale dynamos 

There are some major uncertainties in what exactly are the relevant large scale dynamo 
mechanisms in galaxies and stars. Simulations have enabled us to make close contact between 
simulations and theory. We are therefore beginning to see some significant progress in that 
many of the uncertainties intrinsic to mean field dynamo theory can now be eliminated. 
However, the anticipated agreement between theory and simulations does not yet extend to 
the more realistic cases with strong inhomogeneities and anisotropics, for example where the 
turbulence is driven by convection, supernovae, or by the magneto-rotational instability. 

Crucial to improving the agreement between theory and simulations has been the real- 
ization that the dominant nonlinear feedback comes from the current helicity term. This 
term can produce a "magnetic a effect" (which is aM = ^tJ ■ b/ po in the isotropic case) - 
even if there is no kinetic a effect. The latter is aK = — -Iru; ■ u in the isotropic case, where 
u; = V X M is the vorticity, J = V x 6//io the current density, and po the mean density. 
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Figure 6: Evolution of the energies of the total field (-B^) and of the mean field ), in 
units of -S^q, for runs with non- helical forcing and open or closed boundaries; see the solid and 

dotted lines, respectively. The inset shows a comparison of the ratio (S for nonhelical 

(a = 0) and helical [a > 0) runs. For the nonhelical case the run with closed boundaries is 
also shown (dotted line near {B )/{B'^) w 0.07). Note that saturation of the large scale field 
occurs on a dynamical time scale; the resistive time scale is given on the upper abscissa. 



An example where can be generated - even with ax = - is the shear-current effect of 
Rogachevskii & Kleeorin (2003, 2004). 

Recently, quantitative comparisons between theory and simulations have shown that what- 
ever mean electromotive force {£ = u x b) is produced by the mean field dynamo, this pro- 
duces magnetic helicity in the large scale field. Because of magnetic helicity conservation, 
a corresponding negative contribution in the magnetic helicity of the small scale magnetic 
field must be generated. This results in the production of a small scale current helicity which 
enters in the mean field equations as a magnetic alpha effect. To satisfy magnetic helic- 
ity conservation at all times, the small scale magnetic helicity equation (or cquivalcntly the 
evolution equation for q;m) has to be solved simultaneously with the mean field equations. 
This reproduces quantitatively the resistively slow saturation of helically forced dynamos in a 
periodic box (Field & Blackman 2002, Blackman & Brandenburg 2002, Subramanian 2002). 
The explicitly time-dependent evolution equation for the magnetic a cfi'ect was first derived 
by Kleeorin & Ruzmaikin (1982). Early work by Ruzmaikin (1981) focused attention on 
the possibility of chaotic dynamics introduced by this effect (see also subsequent work by 
Schmalz & Stix 1991 and Covas ct al. 1997), and the first connection between dynamical and 
catastrophic quenching was made by Kleeorin et al. (1995). 

In order to avoid too much repetition with recent reviews on this subject (e.g., Branden- 
burg et al. 2002) we discuss here only in a few words the main consequences of the dynamical 
quenching model. It reproduces quantitatively the resistively slow saturation phase of an 

dynamo in a periodic box (Field & Blackman 2002). Secondly, it reproduces reasonably 
accurately the saturation amplitude and the cycle period of ail dynamos with a sinusoidal 
shear profile in a periodic box. 

In the case of galaxies the dynamical time scale is 10^ yr, which puts rather stringent 
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constraints if one wants to explain microgauss field strengths in very young galaxies that are 
10^ yr old. It is likely that a successful dynamo has to have magnetic helicity fluxes that 
allow the dynamo to get rid of small scale magnetic helicity to allow for a rapid build-up 
of a large scale field that tends to have magnetic helicity of the opposite sign (Blackman & 
Field 2000, Kleeorin et al. 2000). The currently most convincing example where the presence 
of boundaries has been found to be important is in connection with turbulent dynamos that 
work mainly in the presence of shear (Brandenburg et al. 2005). Here, the main effect that is 
thought to be responsible is the so-called shear-current effect where the electromotive force 
has a component proportional to W x J. Here, W is the vorticity of the mean flow and J is 
the mean current density. This effect is technically related to the SI x J effect (e.g. Krause & 
Radler 1980), which is obviously quite distinct from the famous a effect. Another possibility 
is the Vishniac & Cho (2001) mechanism. However, it has not yet been possible to verify 
any of them explicitly. Preliminary mean field calculations suggest, however, that the W x J 
effect produces qualitatively favorable agreement with the direct simulations. 

Figure El shows that in the presence of closed (i.e. perfectly conducting) boundaries, mag- 
netic energy both of the total field and of the mean field saturates at a much lower level than 
with open boundaries {nx B = 0). The inset shows the ratio between the energies contained 
in the mean field to that in the total field. Note that this ratio stays small when the bound- 
aries are closed, but it increases to fairly large values (around 0.7) when the boundaries are 
open. 

4 Conclusions 

In many astrophysical bodies some type of large scale dynamo is likely to operate. This 
dynamo may operate with nonhelical turbulence and shear alone, i.e. without a effect. Small 
scale dynamos, on the other hand, are an extreme case which requires that there is no shear 
and no helicity, for example. Both types of dynamos are vulnerable in their own ways: a large 
scale dynamo requires magnetic and current helicity fluxes in order to be successful, while 
small scale dynamos may require magnetic Prandtl numbers that are not too small. However, 
as we have shown in the present paper, small scale dynamos still work in the compressible 
regime. In fact, it now seems that, once the Mach number exceeds unity, its onset becomes 
independent of the Mach number. This hypothesis has only been tested for small Mach 
numbers, so it would be useful to extent these studies to larger values. However, as we 
have also been able to show, the use of shock-capturing viscosities seems to be a reasonably 
accurate approximation for this purpose. 
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